comparing-maps.Rmdlibrary(Racmacs)Different maps can be compared using the procrustesMap() function where one map is compared to another.
We’ll read in the H3N2 2004 map as an example map
map <- read.acmap(system.file("extdata/h3map2004.ace", package = "Racmacs"))
view(map)We can use the procrustesMap() function to compare the top 3 optimization runs.
pc_run12 <- procrustesMap(
map = map,
comparison_map = map,
optimization_number = 1,
comparison_optimization_number = 2
)
view(pc_run12)pc_run13 <- procrustesMap(
map = map,
comparison_map = map,
optimization_number = 1,
comparison_optimization_number = 3
)
view(pc_run13)